Global 5*0(3) x SO (3) x U(l) symmetry of the Hubbard model on bipartite lattices 
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In this paper the global symmetry of the Hubbard model on a bipartite lattice is found to be 
larger than SO (A). The model is one of the most studied many-particle quantum problems, yet 
except in one dimension it has no exact solution, so that there remain many open questions about 
its properties. Symmetry plays an important role in physics and often can be used to extract useful 
information on unsolved non-perturbative quantum problems. Specifically, here it is found that 
J-H ' for on-site interaction U 7^ the local 5(7(2) x SU(2) x (7(1) gauge symmetry of the Hubbard 

model on a bipartite lattice with N D sites and vanishing transfer integral t = can be lifted to 
a global \SU(2) x SU(2) x U(l)]/Z$ = SO(3) x SO (3) x (7(1) symmetry in the presence of the 
kinetic-energy hopping term of the Hamiltonian with t > 0. (Examples of a bipartite lattice are the 
D-dimensional cubic lattices of lattice constant a and edge length L — N a a for which D = 1, 2, 3, ... 
in the number of sites.) The generator of the new found hidden independent charge global (7(1) 
symmetry, which is not related to the ordinary (7(1) gauge subgroup of electromagnetism, is one 
half the rotated-electron number of singly-occupied sites operator. Although addition of chemical- 
^ , potential and magnetic-field operator terms to the model Hamiltonian lowers its symmetry, such 

terms commute with it. Therefore, its 4 N " energy eigenstates refer to representations of the new 
found global [5(7(2) x 5(7(2) x U(l)]/Z% = SO (3) x SO (3) x (7(1) symmetry. Consistently, we find 
that for the Hubbard model on a bipartite lattice the number of independent representations of the 
group 50(3) x 50(3) x (7(1) equals the Hilbert-space dimension 4 N " . It is confirmed elsewhere 
that the new found symmetry has important physical consequences. 

^ 1 PACS numbers: 02.20.Qs, 71.10.Fd 
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I. INTRODUCTION 

The Hubbard model on a bipartite lattice (for instance one-dimensional, square, cubic, and other D-dimensional 
\^ • cubic lattices) is the simplest realistic toy model for description of the electronic correlation effects in general many- 
I particle problems with short-range interaction. It can be experimentally realized with unprecedented precision in 
7— I ■ systems of ultra-cold fermionic atoms on an optical lattice of variable geometry. For the square and cubic lattices 
one may expect very detailed experimental results over a wide range of parameters to be available^. For instance, 
£\j , recently systems of ultra-cold fermionic atoms describing the Mott-Hubbard insulating phase of the the Hubbard 
model on a cubic lattice were studied 2 . On the one dimensional and square lattices the model has been widely used 
for describing the effects of correlations in several types of materials such as quasi-one-dimensional conductors^ and 
■ high-Tc superconductors^—. 

Unfortunately, most exact results and well-controlled approximations for this model exist only in one dimension 
(ID) 8-10 . Many open questions about its properties remain unsolved. One of the few exact results, which refers to the 
model on any bipartite lattice, is that for on-site interaction U ^ it contains a global SO (4) — [5(7(2) x SU(2)]/Z2 
symmetry. It is associated with a spin 5(7(2) symmetry and a charge ?7-spin 5(7(2) symmetry-ii. We denote the ?7-spin 
(and spin) value of the energy eigenstates by S v (and 5 S ) and the corresponding projection by St = — [N^ — N}/2 
(and 5J = —[N^ — iVj/2). Here denotes the number of lattice sites and N = Nf + iNTj. that of electrons. Our 
notation is particularly appropriate to a _D-dimensional cubic lattice where D = 1, 2, 3, ... for the one-dimensional, 
square, cubic ... lattice, respectively, N a is the number of sites in an edge of length L = N a a, and a is the spacing. 

In this paper we find that for (7^0 the local 5(7(2) x 5(7(2) x (7(1) gauge symmetry of the Hubbard model on 
a bipartite lattice with transfer integral t = <M can be lifted to a global [5(7(2) x 5(7(2) x (7(1)]/Z| = 50(3) x 
50(3) x (7(1) symmetry for the model with t > 0. Indeed, the requirement of commutability with the U/4t ^ 
interacting Hamiltonian replaces the (7 = global 0(4)/Z 2 = [SO (A) x Z 2 ]/Z 2 symmetry by 50(3) x 50(3) x (7(1) = 
[SO (4) x U(l)]/Z 2 rather than 50(4). Here, the factor Z 2 in 50(4) x Z 2 refers to the particle-hole transformation 
on a single spin under which the interacting term is not invariant 12 . In 0(4) /Z 2 and [5(7(2) x 5(7(2) x (7(1)]/Z 2 = 
50(3) x 50(3) x (7(1) the factors 1/Z 2 and l/Zf impose that [S v + S s ] and [S v + S s + 5 C ], respectively, are integers. 
In the latter equations 5 C is the eigenvalue of the generator of the new global (7(1) symmetry found in this paper. 
Our results profit from those of RefJ^ and reveal that such a symmetry becomes explicit, provided that one describes 
the problem in terms of rotated electrons. Those are generated by any of the unitary transformations considered in 
that reference, which refer to U/4t > values and can be trivially extended to U/4t ^ values. The global symmetry 



2 



found here refers to the latter U/At range. 

The paper is organized as follows. The model and the electron - rotated-electron unitary transformations are the 
subjects of Section II. In Section III a global 5*0(3) x 5*0(3) x U(l) symmetry is established for the Hubbard model 
on a bipartite lattice with U/At ^ 0. Finally, Section IV contains the concluding remarks. 

II. THE MODEL AND A SET USEFUL ELECTRON - ROTATED-ELECTRON UNITARY 

TRANSFORMATIONS 

On a bipartite lattice with spacing a, = [N a ] D sites, N a even, N a /2 odd, L — N a a, and spatial dimension 
D < N a the Hubbard model is given by, 

U 

H = f+-{N°-Q}; f=-tJ2 E [4,* C V- + h - c \ ; Q = E E h n° (1 - • (1) 

(Fjry) <T=f4 j=l<r=t4 

Here T is the kinetic-energy operator with first-neighbor transfer integral t, which can be expressed in terms of the 
operators, 

To = ~ E E^i - CT 4j,a^'^%'.-^ + ( 1_ % -^)4 3 -,a^'^( 1- ^' ,-ff) + ^- C -]' 

2+1 = - E E C ?j,a C r,',cr (1 ~ %y,-g) , 

(f.,-fy) o- 

^-i = - E E^ 1 - CT )4 ) <7 c v^" f 'i' -°-> ( 2 ) 

as T = t [Tq + T+i + T_i]. While the operator To does not change electron double occupancy, the operators T + i and 
f—i do it by +1 and —1, respectively. In the above equations n?. a — ct ^c?^, ±<j refer to opposite spin projections, 

and the operator Q counts the number of electron singly occupied sites. Hence the operators, 

D = \[N-Q]; D h = 7^[N h - Q] ; Q t = \ [Q + (iV t - jfy] ; Q x = \ [Q - (jV t - Nj] , (3) 

count the number of electron doubly occupied sites, unoccupied sites, and spin a =f, \ singly occupied sites, respec- 
tively. Moreover, TV = ft a and N a = J^fli n?,,* where N h = [2iVf - N] , ft* = [iVf - ft±] , and ft£ = [N° - N t ] . 

For simplicity let us consider that U /At > and let {jvl/oo)} be a complete set of 4^ energy eigenstates for U /At — > 
co. There is exactly one unitary operator V = V(U/At) such that for the value of U /At > under consideration each 
of the A N " states \^u/it) = V^l^oo) is generated from the electronic vacuum by the same occupancy configurations 
of rotated electrons of creation operator ct as the corresponding U/At — > oo energy eigenstate in terms of electrons. 
The rotated-electron creation and annihilation operators are given by, 

4„ CT = ^ t 4,a^; C-r 3 ,a = V^ C F] , a V- h? j>a = ct ^ C^., CT . (4) 

Rotated-electron single and double occupancy are good quantum numbers for U/At > whereas for electrons such 
occupancies become good quantum numbers only for U/At oo. Therefore, V = V(U/At) becomes the unit operator 
in that limit. The unitary transformation associated with the operator V is of the type studied in Ref.— . There is one 
of such transformations for each choice of U/At — > oo energy eigenstates. Similar results are obtained for U/At < 0. 

We introduce the operator O — OV. It has the same expression in terms of rotated-electron creation and 
annihilation operators as O in terms of electron creation and annihilation operators. Here V — V . Note that within 
our representation both the notations referring to marks placed over letters being a caret W or a tilde L denote 
operators. Such notations are useful for operators for which W — L such as the general operators and 0. Indeed, 
then they imply the equivalent relations = V* O V and = VOV^. (Here we have used that V — V.) When 
O O our convention is that in general the expression of the operator O in terms of electron creation and annihilation 
operators is simpler than that of = V* V in terms such operators. This then implies that the expression of 
in terms of rotated-electron creation and annihilation operators is simpler than that of O = V in terms of the 
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same rotated-electron operators. (An exception are the electron operators of Eq. (gjl, which denote by ct a and cp^ a 
rather than by ct. a and c? j}(7 , respectively.) 
Any operator O can be written as, 

6 = vov^ = o + [6, s} + ~[[6, s], §] + ... = vov* = o + [6, §} + ^[[6, s], §} + ..., (5) 

where = e s , V — e~ s , and S = S. That S and V have the same expression both in terms of electron and 
rotated-electron creation and annihilation operators justifies that O = V OV^ = V O in Eq. |5j). Importantly, it 
follows from the results of Refj^ that for each electron - rotated-electron unitary transformation and corresponding 
unitary operator V of the type considered in that reference the operator S has a uniquely defined expression. For 
any of such transformations that unknown expression of S involves only the kinetic operators To, 7+ 1, and T_i of Eq. 
(|2|) and numerical U/At dependent coefficients. For U/At ^ it can be expanded in a series of t/U. Importantly, the 
corresponding first-order term has a universal form for all electron - rotated-electron unitary transformations of the 
above-mentioned type, which readsi^, 



+ o(t 2 /u 2 ) = s = -^ T +1 -r_, +c7(/-/r-i. ((» i 



(The form of our relation = e s justifies the extra minus sign in the S and S expressions given here, relative to 
those of ReiM.) 

Furthermore, for any unitary operator V of the above type, —S can be written as —S — S(oo) + AS. Here S(oo) 
corresponds to the operator S(l) for I = oo defined in Eq. (61) of Ref.— and AS has the general form provided in Eq. 
(64) of that reference. For each specific transformation and corresponding choice of U/At — > oo energy eigenstates there 
is exactly one choice for the numbers D^(m) in that equation, (k = 1, 2, ... refers to the number of rotated-electron 
doubly occupied sites.) 

Since V is unitary, the operators ct a and cp ja have the same anticommutation relations as ct ct and c? j(J . The 



a electron number operator N a = "£^=1 ^f-,a- equals the corresponding a rotated-electron number operator N a = 



^2^=1 npj ,a ■ As a result, it remains invariant under V, so that [iV^V] = [N a ,S] = 0. (See equation (O such that 
[N a , S}=0 for 6 = N a and 6 = N a .) 



III. THE GLOBAL SO (3) x SO (3) x 1/(1) SYMMETRY FOR U/At / 

A. Global symmetry of the Hubbard model on a general bipartite lattice 

The local SU(2) x SU(2) x U(l) gauge symmetry of the Hamiltonian (TTJ) for U/At —> ±oo considered in Ref^ 
becomes for finite \U/At\ > values a group of permissible unitary transformations. It is such that the corresponding 
local U(l) canonical transformation is not the ordinary U(l) gauge subgroup of electromagnetism. Instead it is a 
"nonlinear" transformation 12 . Following the unitary character of V — V, one can either consider that, 

H=VHV^ = VHV^ = H+[H, §] + i[[H, S], §] + ..., (7) 

is the Hubbard model written in terms of rotated-electron operators or another Hamiltonian with an involved expres- 
sion and whose operators ct.^ and c? j(7 refer to electrons. According to RefJ^, the latter rotated Hamiltonian is built 

up by use of the conservation of singly occupancy 2S C = (Q) by eliminating terms in the t > Hubbard Hamiltonian. 
That is done so that S c is an eigenvalue of the following one-half rotated-electron singly-occupancy number operator 
associated with the operator S c = Q/2, 

Sc^^QV= l -Q= l -Y^Y. ^(l-fifw). (8) 

3 = 1 <r=t4 

Here h? jt(T — h? jM V — ct , a c? ja is the operator given in Eq. Q. According to the studies of RefJ^, this can be 

done to all orders of t/U provided that U/At ^ 0. In the context of Ref. 1 -, this is equivalent to compute rotated 
"quasicharge" fermions whose number exactly equals [N^ — 2S C ]. 
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The "rotated" Hamiltonian H = V* H V commutes with the six generators of the SO{A) symmetry. Thus the 
Hubbard model H commutes with both such generators and corresponding six other operators with the same expres- 
sions when written in terms of rotated-electron operators. Consistently with Eq. ([5]), this just means that the six 
generators of the f?-spin and spin algebras commute with V . To reach this result we have profited from the expression 
of the operator S only involving the three kinetic operators given in Eq. ((2]). We have then calculated the following 
commutators, 



[&a,Tl] 



[S a ,Ti] 



0: 



rj, s , 2 = 0, ±1 



(9) 



Although the algebra involved in their derivation is cumbersome, it is straightforward. Therefore, we omit here the 
corresponding details. The vanishing of the commutators ([9]) implies that the six generators of the T^-spin and spin 
algebras commute with V, 



[S z a ,V] = [SlV] = [S a ,V] = 0; 



77, s. 



(10) 



This confirms that for such six operators all operator terms on the right-hand side of Eq. ([5]) containing commutators 
vanish so that O = O for O being any of such operators. Hence they have the same expression in terms of electron 
and rotated-electron operators and read, 





i 


[N?-N} = -~[N?-N]; S* s = - 


~[^t--^] = - 








N a 




N a 




= E 




S v = J2 er l *- ? i 


Cf 3 t °rji — E e 




i=i 


3=1 


3 = 1 


3=1 












si 


= E 


4ji c ^t = E %i c? jt i Ss = E 


4 it ^=E4 3 


t C 3\ I ' 




3=1 


3=1 3=1 


3 = 1 





(11) 



where the vector tt has Cartesian components 7? = [tt, tt, ...]. For instance, for the model on the ID, square, and cubic 
lattices those read 7T, [7T, tt], and [7r,7r,7r], respectively. 

In addition, we have evaluated the commutators of the three components of the momentum operator P with the 
three operators of Eq. ([2]). Again all such commutators vanish, so that the momentum operator commutes with V . 
Use of Eq. <j5j> then implies that such an operator reads, 



E 

o=t,4 % 



E H 



E E 

<T=t,i k 



(12) 



Again all operator terms on the right-hand side of Eq. containing commutators vanish for O being any of the 

three operator components of P, so that P = P. 

According to the studies of RefJ^, the SU(2) x SU(2) x U(l) Lie group and its local generators can be represented 
by the 4x4 on-site matrix x?. provided in Eq. (7) of that reference and matrices op. appropriate to these generators. 

Their entries are given through polynomials of electron operators of the general form Xf>, = v Xf.jji rhp-i'j = 

Tr {xp rhf.) and Op. = v op. i i> 'fhr-,1',1 = ^(ofj rhf), respectively. Here the operator matrix rh?. has the same 

form as the operator matrix fhp j = V' rhp j V , but with the rotated-electron operators replaced by electron operators. 
The operator matrix fhp. plays an important role in our studies. It reads, 



1 - n?j ,t - n?i 4 + np } ,t n? 3 4 cp. 4 c ?j 4 (1 - n ?j 4) cp. 4 (1 - np. 4) cp j 4 



£1^(1-^.4) 

%4( 1_ "^'t) 



-Tiff. 



^ ,t 4- 4 >t ( 1 " ,l) 4, 4 ~ c ?o 4 



™r 3 4 cr, 4 cl ^ c ?j 4 fifj 4 (1 - rift. 4) 



(13) 



As described in RefJ^ for the polynomial Op. , one can as well introduce a general polynomial operator Op. of rotated- 
electron operators of the general form, 



(14) 



U' 
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Lifting the local ?7-spin and spin SU(2) x SU(2) gauge symmetry of the Hubbard model on a bipartite lattice for 
U/4t = ±00 to a global [SU(2) x SU(2)]/Z 2 = 50(4) symmetry of that model for U/At ^ is simply accomplished by 
summing over the sites the six local generators Op of the SU(2) x SU(2) sub-group of the SU (2) x SU(2) x U (1) 
Lie group. It follows from the equalities of Eq. that the six generators given in that equation can be represented 

by polynomials of electron and rotated-electron operators of the same form, Ylj=i Ofy = 2j=i Op j ■ This holds in 

spite of except for U/At — > ±00 the corresponding local generators Op, and Op. being different operators, Op j ^ Op y 

Indeed, the local generators Op. do not in general commute with the unitary operator V. This follows from rhp^ y \ ^ 

rhpjj'j, where rhp.^^ and mp jf v t i appear in the expressions Op, = Y,i,i> °r,,i,i> mp jt i> t i = Tr (op, m^.) and ((HI) of Op d , 
respectively. However, the matrix op. appearing in these two expressions is the same. For the six local generators 
associated with the generators (ITU of the global 5*0(4) symmetry it reads, 



-1/2 

1/2 









-1/2 
1/2 



(15) 



for the ?7-spin and spin diagonal generators and 
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- 
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1 


0. 



(16) 



plus their two hermitian conjugates for the r^-spin and spin off-diagonal generators. 

Now for the "rotated" Hamiltonian H = V^HV = V^HV& local SU(2) x SU{2) x {7(1) gauge symmetry occurs 
for U/At — > ±00 as well. Alike the original Hamiltonian, H has a global SO(A) symmetry whose generators are 
obtained as above. In addition, a similar procedure can be used to lift the local U(l) gauge symmetry to a global 
symmetry of the "rotated" Hamiltonian for t > and U/At ^ 0. Indeed, through the polynomial of rotated-electron 
operators given in Eq. (fT4")) . the local generator of the "nonlinear" local U(\) gauge symmetry can be represented by 
a 4 x 4 matrix given by, 







1/2 

1/2 



(17) 



This local generator refers to rotated-electron single occupancy 2S C . The use of such a matrix op. in the polynomial 

Opj of Eq. ([TJ} leads for U /At ^ to a sum of polynomials J2j=i Opj ■ It exactly equals expression © of the generator 
of the global U(l) symmetry whose eigenvalue S c is one half the number of rotated-electron singly occupied sites 2S C . 
The trivially related operator S% = [D + D h }/2 of eigenvalue = [iVf/2 - S c ] can also generate such a global 
symmetry of the "rotated" Hamiltonian H = V' H V = H V. When written in terms of local polynomials as 
"52^=1 Opj , its corresponding matrix op j reads, 



Op 



1/2 

1/2 







(18) 



This generator refers to rotated-electron double occupancy plus non occupancy [2^^ — 2S C ). However, pA^f — 25 c ] 
and rotated-electron single occupancy 2S C are not independent. Hence the operators S c = Q/2 of Eq. associated 
with the matrix opj of Eq. (fTTj) and S*' 1 = [D + D h }/2 associated with the matrix op j of Eq. (TI5)) refer to two alternative 
representations of the generator of the global 17(1) symmetry of the "rotated" Hamiltonian under consideration. 

The main point is that a global U(l) symmetry in the "rotated" Hamiltonian H — V* H V — V' H V for t > 
and U/At ^ must also be a global U(l) symmetry, which is hidden in the original model H = V HV> = V HV^ 
of Eq. (|7J| . Indeed, for the latter original model the generator flSJ) refers to one half the number of rotated electrons 
rather than electrons. And in contrast to the six generators (|TT|) of the global SO (A) symmetry, the number of rotated 
electrons operator does not commute with the unitary operator V. Other related operators D, D h , and Q ai which for 
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U/At 7^ also commute with the Hamiltonian (JXJ) yet do not commute with V, are obtained by rotating the number 
operators D, D h , and Q a , respectively, given in Eq. ([3]). For the A^-site problem only for rotated electrons does 
single and double occupancy remain good quantum numbers for finite \U/At\ > 0, whereas for electrons single and 
double occupancy are conserved only for \U/At\ — > oo. This is why the generator (jHJ does not commute with V. 

Since is even, both [S n + S%] and [S* + Sf ] are integers. Their relation to 5 C is such that 25 c and [N^ — 25 c ] give 
the number of spin-1/2 spins of the rotated electrons that singly occupy sites and the number of r?-spin-l/2 "77-spins", 
respectively. The former number equals that of rotated-electron singly occupied sites and the latter number that of 
rotated-electron doubly occupied sites (down ry-spins) plus rotated-electron unoccupied sites (up ry-spins) , respectively. 
Therefore, [SJ + S% + S c ] must also be an integer. This justifies why for U /At ^ the global symmetry of the model 
Q on a bipartite lattice is that of the group [SU{2) x SU{2) x U(l)]/Z$ = 50(3) x 50(3) x £7(1) rather than 
SU(2) x SU(2) x U(l). The global £7(1) symmetry remained hidden because in contrast to the six generators 

one has that the generator S c = YLj=i 0?j is except for \U/At\ — > 00 different from the operator Ylj=i Or } ■ (Note that 
the matrix op. is given by Eq. (|17[) both in the Op. and Op. expressions.) Indeed, when written in terms of electron 
creation and annihilation operators the expression of the generator 5 C is for \U/At\ finite involved, consisting of an 
infinite number of operator terms, 

S c = £ W 6 P] V = J2 (Op, + , ^ ] + - [[6 fj , S+ ], St ] + ...) , (19) 
rather than merely by Ylj=i Opy Only for \U/At\ — > 00 one has that the commutator [Op j , S'\ = vanishes, so that 

Sc = S c = ^2j=i ■ 



B. Consistency between the global symmetry and the Hilbert space dimension 

Addition of chemical-potential and magnetic-field operator terms to the Hamiltonian (JTJ) lowers its symmetry. 
However, such terms commute with it. Therefore, the global symmetry being [SU(2) x SU(2) x C/(l)]/Z| = 5*0(3) x 
50(3) x U(l) implies that the set of independent rotated-electron occupancy configurations that generate the model 
energy eigenstates generate state representations of that global symmetry for all values of the electronic density 
n and spin density m. It then follows that the total number of such independent representations must equal the 
Hilbert-space dimension 4 . Here we show that for the Hubbard model on a bipartite lattice the number of 
independent representations of the group 50(3) x 50(3) x U(l) is indeed 4 . In contrast, the number of independent 
representations of the group 50(4) is for that model found to be smaller than its Hilbert-space dimension A N " . This 
is then consistent with the global symmetry of the Hubbard model on a bipartite lattice being larger than 50(4) and 
given by 50(3) x 50(3) x £7(1). 

The rotated-electron occupancy configurations involving the (i) singly occupied and (ii) unoccupied and doubly- 
occupied sites are independent. They refer to the state representations of the spin SU{2) symmetry M 8 = 25 c spin-1/2 
spins and 77-spin SU{2) symmetry M v — 25^ ?7-spin-l/2 ry-spins, respectively. Indeed, concerning the 77-spin SU{2) 
representations the rotated-electron doubly occupied sites and unoccupied sites play the role of down and up 7y-spin-l/2 
?7-spins, respectively. In turn, the U(l) symmetry state representations refer to the relative occupancy configurations 
of the 25 c rotated-electron singly-occupied sites and 25^ rotated-electron unoccupied and doubly-occupied sites. For 
U/At ^ 0, the Hilbert space can then be divided into a set of subspaces with fixed S v , S s , and 5 C values and thus with 
the same values M v = 25^ of 77-spins and M s = 2S C of spins. The number of SU(2) x SU{2) state representations 
with both fixed values of S v and 5 S , which one can generate from M v ?7-spin-l/2 77-spins and M s spin-1/2 spins, reads 
M(S v ,Mi).M{S. t M.). Here, 

Mil., M„, = p*. + 1) { ( M Jl J - ( MJ2 \ _,)}, (20) 

for a — 77, s. If for U/At ^ the global symmetry of the model was SO(4), then the dimension of such a subspace 
would be A/"(5 T? , M v ).J\f(S s , M s ) and the sum of all sub-space dimensions would give the Hilbert-space dimension 
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4 Na . However, we confirm below that such a sum obeys the inequality, 

EEEI1 2 l -ms a ,M a ) 

M»=0 S„=0 S.,=0 a=rj,s 
JVf M,/2 M./2 i->[2S 0( -(M,-A!f J) )/2]l 

= E E E II [ ( } 2 W(5 a ,M a )<4<\ (21) 

and thus corresponds to a dimension smaller than . Note that M n = [N® — M s ] so that the numbers M v and 
M 8 are not independent. Therefore, the sums over M s = 0, 1, and M v — 0, 1, are indeed alternative, as 
given in Eq. ([2"TT) . 

In turn, that the model global symmetry is larger than ,50(4) = [SU(2) x SU(2)]/Z 2 and given by 50(3) x 50(3) x 
U(l) = [SU (2) x SU(2) x C/(l)]/Z| requires that one accounts for the number of representation states of the extra 
global C/(l) symmetry. For U /At ^ it has in the subpaces considered here, 



representation states. Thus rather than M{S v ,M ri ).M{S s ,M s ) each of such subspaces has a larger dimension, 

d{S n ,S B ,S c ) = d c M{S v ,N? -2S C )M{S.,2S C ). (23) 
By performing the sum over all subspaces, one then finds indeed in Appendix A that, 



[JVf /2] [JV a D /2-S c ] 5 C 

= E E E 

s c =o s„=o s 3 =o 
x (25, + 1) {2S S + 1) 
N? 



(N a D 
\2S C 



a.—r}.s 
D 



ND/2 - S v + 5, 



A^/2 - 5, - 5 S 



1 



TVf/2 - 5„ + 5 



s„=o s,=o 

N?/2 -S v - Ss) \N°/2 - S v - S 8 - 2 
) + ( 

+ \N?/2 -S v + S s -1 



(24) 



Finally, except that a factor of one in each term of the two alternative sums of Eq. (f2lj) is replaced by the 
dimension d c in the sum of Eq. (|24p , such sums are identical. This is confirmed by transforming the sums over M s = 
0, 1, and M v = 0, 1, N° of Eq. ((HJ) in sums over 5 C = M s /2 = 0, 1/2, /2 and 5 C = N» /2 - M n /2 = 
iVf /2, iV a D /2 - 1/2, 1/2, 0, respectively, and accounting for that 25 c can be written as 25 c = N^/2- (M r) - M s ) /2 
where /2 is an odd integer number. That the dimension of Eq. (|22|) obeys the inequality d c > 1 then implies the 
validity of the inequality given in Eq. (|2"T1) . 



C. Relation of the global symmetry to the exact solution of the ID model 

In the particular case of the bipartite ID lattice the Hubbard model has an exact solution^—. Since the global 
50(3) x 50(3) x U(l) symmetry found here refers to ID as well, it must be related to that exact solution. Such 
a solution refers to the ID Hubbard model in the subspace spanned by the highest-weight states (HWSs) or lowest- 
weight states (LWSs) of both the 77-spin SU(2) and spin SU(2) algebras. The model energy eigenstates that are HWSs 
or LWSs of these algebras are often called Bethe states. In order to clarify such a relation, rather than the so called 
coordinate Bethe ansata&S, it is convenient to consider the exact solution of the problem by the algebraic operator 
formulation of RefJ^. Within the latter formulation the HWSs or LWSs of the ?7-spin and spin algebras are built up 
in terms of linear combination of products of several types of annihilation or creation fields acting onto the hole or 
electronic vacuum, respectively. 

The algebraic formulation of the Bethe states refers to the transfer matrix of the classical coupled spin model, 
which is the "covering" ID Hubbard models. Indeed, within the inverse scattering metho d 10 ' 16 the central object 
to be diagonalized is the quantum transfer matrix rather than the underlying ID Hubbard model. The transfer- 
matrix eigenvalues provide the spectrum of a set of conserved charges. The diagonalization of the charge degrees of 
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freedom involves a transfer matrix associated with a charge monodromy matrix of the form provided in Eq. (21) of 
RefiiS. Its off-diagonal entries are some of the creation and annihilation fields. The commutation relations of such 
important operators are given in Eqs. (25), (40)-(42), (B.1)-(B.3), (B.7)-(B.ll), and (B.19)-(B.22) of that reference. 
The solution of the spin degrees of freedom involves the diagonalization of the auxiliary transfer matrix associated 
with the spin monodromy matrix provided in Eq. (95) of Ref. 10 . Again, the off-diagonal entries of that matrix play 
the role of creation and annihilation fields, whose commutation relations are given in Eq. (98) of that reference. The 
latter commutation relations correspond to the usual Faddeev-Zamolodchikov algebra associated with the traditional 
ABCD form of the elements of the monodromy matrix^. It also applies to the ID isotropic Heinsenberg model, whose 
global symmetry is SU(2). Consistently, at half filling and for large U /At values the latter model describes the spin 
degrees of freedom of the ID Hubbard model. In turn, the above relations associated with the charge monodromy 
matrix refer to a different algebra. The corresponding form of that matrix is called ABCDF by the authors of RefJ^. 

The main reason why the solution of the problem by the algebraic inverse scattering method^ was achieved only 
thirty years after that of the coordinate Bethe ansatz^S is that it was expected that the charge and spin monodromy 
matrices had the same traditional ABCD form, found previously for the related ID isotropic Heinsenberg models. 
Indeed, such an expectation was that consistent with the occurrence of a spin SU(2) symmetry and a charge (and 
77-spin) SU(2) symmetry known long agoii, associated with a global 50(4) = [SU{2) x SU(2)]/Z 2 symmetry. If that 
was the whole global symmetry of the ID Hubbard model, the charge and spin sectors would be associated with the 
77-spin SU{2) symmetry and spin SU{2) symmetry, respectively. A global 5*0(4) = [SU{2) x SU(2)]/Z 2 symmetry 
would then imply that the charge and spin monodromy matrices had indeed the same Faddeev-Zamolodchikov ABCD 
form. 

However, all tentative schemes using charge and spin monodromy matrices of the same ABCD form failed to achieve 
the Bethe-ansatz equations obtained by means of the coordinate Bethe ansata&S. Fortunately, the problem was solved 
by Martins and Ramos, who used an appropriate representation of the charge and spin monodromy matrices, which 
allows for possible hidden symmetries. Indeed, the structure of the charge and spin monodromy matrices introduced 
by these authors is able to distinguish creation and annihilation fields as well as possible hidden symmetries. 

Our results refer to the Hubbard model on any bipartite lattice. Hence for the particular case of the bipartite 
ID lattice they show that the hidden symmetry beyond 50(4) is the charge global U(\) symmetry found in this 
paper. Our studies reveal that for U /At > the model charge and spin degrees of freedom are associated with U(2) = 
SU(2) x U(l) and SU(2) symmetries, rather than with two SU{2) symmetries, respectively. The occurrence of such 
charge U{2) — SU{2) x U(l) symmetry and spin SU(2) symmetry is behind the different ABCDF and ABCD forms 
of the charge and spin monodromy matrices of Eqs. (21) and (95) of Ref^, respectively. Indeed, the former matrix 
is larger than the latter and involves more fields than expected from the model global 50(4) = \SU(2) x SU(2)]/Z2 
symmetry alone. This follows from the global symmetry of the model on the ID and other bipartite lattices being 
50(3) x 50(3) x U(l) = [SU(2) x C/(2)]/Z| rather than 50(4) = [SU(2) x SU(2)]/Z 2 , as found in this paper. Hence 
our general results for the Hubbard model on a bipartite lattice are consistent with the algebraic operator formulation 
of its exact solution for the particular case of the ID lattice^. 

IV. CONCLUDING REMARKS 

On a square lattice, the Hubbard model is one of the most studied condensed-matter quantum problems. Further- 
more, on any bipartite lattice it is the simplest realistic toy model for description of the electronic correlation effects 
in general many-electron problems with short-range interaction. Therefore, that the global symmetry of the Hubbard 
model on a bipartite lattice is larger than 50(4) and given by 50(3) x 50(3) x U(l) is an important exact result in 
its own right. Furthermore, the new found global symmetry is expected to have important physical consequences. 

The studies of Refp^ 7 - on the Hubbard model on the square lattice use a description in terms of quantum objects 
related to the rotated electrons. The introduction of such a description involves the global symmetry found in this 
paper and corresponding transformation laws under a suitable electron - rotated-electron unitary transformation of 
the type considered here and in RefJ^. The spinless c fermion, spin-I/2 spinon, and ?7-spin-I/2 r7-spinon operators of 
such a description are a generalization to U /At > of the U /At 3> 1 "quasicharge" , spin, and "pseudospin" operators 
of Ref.™, respectively. The former quantum objects emerge from a suitable electron - rotated-electron unitary trans- 
formation. Their operators have the same expressions in terms of rotated-electron creation and annihilation operators 
as those of Ref in terms of electron creation and annihilation operators, respectively. The occupancy configura- 
tions of the spinless c fcrmions, spin-1/2 spinons, and ?7-spin-l/2 77-spinons generate a set of complete states. Those 
correspond to representations of the U(l), spin SU(2), and 77-spin SU(2) symmetries, respectively, associated with 
the three dimensions of Eq. (|23[) and the global symmetry found in this paper. 

The square-lattice quantum liquid introduced in Ref. 17 contains the one- and two-electron excitations of the Hubbard 
model on a square lattice. At hole concentration x = [N% - N]/N% = 0, U/At « 1.525, and t « 295 meV it is found 
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in that reference to quantitatively describing the spin- wave spectrum observed in the parent compound L^CuO^i^. 
A system of weakly coupled planes, each described by the square-lattice quantum liquid of RefJ^, is the simplest 
realistic toy model for the description of the role of correlations effects in the unusual properties of the cuprate hight- 
temperature superconductors^—. After addition of such a weak three-dimensional uniaxial anisotropy perturbation, 
the Hamiltonian terms that describe the fluctuations of two important pairing phases are for intermediate U/4t 
values found to have the same general form as the microscopic Hamiltonian given in Eq. (1) of Ref. 7 . The main 
difference is that the electron creation and annihilation operators appear replaced by rotated-electron creation and 
annihilation operators, respectively. Evidence is provided elsewhere that such a quantum liquid has for a well-defined 
hole-concentration range a long-range superconducting order. In addition, it seems indeed to contain some of the 
microscopic mechanisms behind the unusual properties of the hole-doped cuprate hight-temperature superconductors. 
It is commonly understood that Hamiltonian symmetries by themselves are not sufficient to prove that a particular 
symmetry is broken in the ground state. However, the symmetry of the action that describes the fluctuations of 
the phases of such a quantum liquid and of that of Refi£ is a global superconducting U(l) symmetry. In the case 
of the former quantum liquid the representations of such a (7(1) symmetry are generated by c fermion occupancy 
configurations. Thus it is directly related to the original model hidden global U(l) symmetry found in this paper, 
whose representations are also generated by c fermion occupancy configurations. Such a preliminary result seems to 
confirm the important role plaid by the hidden U(l) symmetry of the global 5*0(3) x 5*0(3) x U(l) symmetry found 
in this paper for the Hubbard model on a bipartite lattice. 



We thank Alejandro Muramatsu, Tiago C. Ribeiro, and Pedro D. Sacramento for illuminating discussions, Tobias 
Stauber for calling our attention to Rcf.— and for discussions, and the ESF Science programme INSTANS and the 
Portuguese PTDC/FIS/64926/2006 grant for support. 



In this Appendix we perform the subspace-dimension summation of Eq. (|24[) that runs over 5 C , 5^, and S s integer 
and half-odd- integer values. For simplicity here we consider the square lattice so that D = 2 in Eq. (|24[). yet the 
derivation proceeds in a similar way for any other D-dimensional cubic lattice where D — 1,2, 3,.... More generally, 
the sum-rule (|2"4"]l is valid for the Hubbard model on any bipartite lattice. The subspace dimensions have the form 
dr ■ ri Q =^ s A/"(5 Q ,M a ) given in Eq. (|23|) where JV(S a , M a ) and d r are provided in Eqs ((20)) and (|22|) . respectively. 
Recounting the terms of Eq. (|24p . one may choose to be the independent summation variable what gives, 
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Appendix A: Subspace-dimension summation 



Nl/2[Nl/2-S c \ S a 




s c =o 



5„=0 S s =0 




(Al) 



1 + (-1) 2 ( S V+ S *) 1 + (_1)2(S S +S C ) 



2 2 



s„=o s s =o s c =s s 



One can then rewrite the summation C|24[) in the form, 



Nl/2 [N%/2-S n ] 



1 + (-1)2(S,+S 3 ) 



2 



(25^ + l)(25 s + l) x E(Sr,,S.), 



(A2) 



s„=o s s =o 
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where E{S V ,S S ) denotes the S^ and S s dependent summation over S c as follows, 

N a 



£{S V ,S S ) = E 



" 1 + (_i)2(s s +s c ) r N 2 



s c =s s 



Nl - 2S r . 



N 



2S C ; 

Nl - 2S r 



-f- S c Srj 1 



25 c 



2 £" i+(-i) a < a - +g -> ^. 



s c =s s 



1 



25 c 
/Sc <Sg 1 

1 



(S c - S s )\ {S c + S s )\ (S c -S s - 1)! (S c + S s + 1)! 

1 



(JV2/2 - 5 C - 5,)! (iV Q 2 /2 - S c + 5,)! (iVf/2 - S c - S v - 1))! {Nl/2 -S c + S v + 1)! 



(A3) 



In order to evaluate £{S V , S s ) it is useful to replace the variable S c by fc = S c — S s . To simplify the notation we 
then introduce, 



S — Sjj + S s — S(Sjj, S s ) ', 2? — Srj — S s — "D(Sjj, S s ) ■ 
Due to the parity factor, in the summation over fc only the terms with k integer survive so that, 



(A4) 



z = E W 



k=0 



1 



1 



k\{s-v + k)\ (k-iy.(s -v + k + iy. 
i 



i 



E N > 



(JV2/2 - S - fc)! (iV|/2 + 2) - k)\ {Nl/2 -S-k-l)\ {Nl/2 + V-k + l)\ 
1 1 



k=0 



M(S-V + k)\ {Nl/2 -S-k)\ {Nl/2 + X> - k)\ 
1 1 



fc! (5 - V + k)\ {Nil 2 - S - k - 1)! {Nl/2 + V - k + 1)! 
1 1 



(fc - 1)! {S - V + k + 1)! {Nl/2 - S - fc)! (JV2/2 +£> - fc)! 
1 1 



{k-l)\{S -V + k+ 1)1 {Nl/2 -S-k -1)1 {Nl/2 + V - k + 1)! 

where now the fc summation runs over integers only. 

In order to perform the summation (|A5[) we rearrange the terms as follows, 



(A5) 



z= E 



i 



k=Q 



N: 



l/2 + V\f Nl/2 



V 



Nl/2 + V + l\f Nl/2 



fc 



-V - 1 



{Nl/2 + V)\{Nl/2-V)\ 
1 

{Nl/2 + V +1)1 {Nl/2- 25-1)! 
1 

(7V2/2 + 2? - 1)! (7V2/2 - V + 1)! V fc-1 J\N%/2-S-k 



Nl/2 + V\f Nl/2-V 



k J\Nl/2-S-kJ V )\ N l/ 2 -S -k-1 



\Nl/2 -S-k-1 



Nl/2 + V-l\ (Nl/2 -V + l 



>Nl\, 



(A6) 



or 



Ul/2 - V 



E 

fe=0 



JV 2 /2 + P\/ N 2 J2-V 
k ) \Nl/2 -S-k 



Nl/2 + V\f Nl/2 - V 
k-l J\N£/2-S-k-l 



Nl j y 

N^/2-V-l) ^ 



iV 2 

,7V 2 /2 - V + 
Next, by using the identity, 



fc=0 



fc=0 



Nl/2 + V+l^f Nl/2-D J 



E 

fc=0 



k \N -k 



N%/2 


-5 


- fc 






+ r 


Nl/2 


-5 






•A + 




)■( 


TV 


1 



we carry out separately the summations in expression (|A7j) . what gives 



fT 5 /iV Q 2 /2 + P\/ iV 2 /2 



\Nl/2-S-k) Wa/2-S 



E 

fc=0 



Nl/2 + V 
k-l 



Nl/2 - V 
N*/2-S-k-l 



Nl/2-S-l 

E 

fe=i 



jV 2 /2 + 2?\/ N 2 J2-V 
k-l )\N^/2-S-k-l 



/ Nl/2 + V\ Nl/2-V 



E 

fc'=0 

NI/2-S-2 



\N%/2-S -2-k' 



and 



2^t \ k I V 7V2 /9 - 

fc=0 v 



2/0_ 2? _ 1 

V iV 2 /2-S-fc-l 



/ Nl/2 + V + l\( 7V 2 /2-2?-l 



E 

k=0 

Nl 

m/2-s-i 



\Nl/2-S-l 



N 2 /2-S 

E 

k=0 



Nl/2 + V-l\ fNl/2-V + l 
k-l J \n*/2 -S-k 



jy2 y2 5 

/ 7V 2 /2 + P-l\ fNl/2-V + l 
k-l ) \N*/2 -S-k 



E 

fc=i 

iV 2 /2-<S-l 



fc'=0 

iV 2 

v iV 2 /2 - 5 - 1 

Introducing these results in expression (|A7|) for i7 leads to, 

TV 2 \ [/ iV 2 \ / Nl 
y N%/2 - V) [Ua 2 /2 - S) + Ua 2 /2 -S-2 
Nl \ \( N 2 
,Nl/2 -S-l 
E(5,D). 



v /iV 2 /2 + 2?-l\ / iV 2 /2-2? + l 

\ k' )\Nl/2-S-l-k' 



EySrf) S s ) 



NZ/2-V + 1J \Nl/2-V -1 
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Expression (|A2[) for Aftot can now be rewritten as, 



1 + (-l)2(^+S s ) 



[Af 2 /2-S„] 

Mo, = E E 

s„=o s s =o 

,JV2/2 - (S„ - S s ) 
x N%/2 - (S^+ S s ) - 1 



(25 n + l)(2S' s + l) x 



Nl/2 - (S v + S s )) + \Nl/2 - (S v + S a ) - 2 

Nl/2 - (S v - S s ) + l) + U a 2 /2 - (S v -S s )-l 



(A14) 



where the summations run over both integers and half-odd integers. The use of the notation (|A4[) then allows rewriting 
(IA14[) in compact form, 



Nl/2 [Nl/2-S v ] 



Mo* = E E 



i + (-i) 



2S 



-(S + V+ l)(S-V+ 1) xS(5,P) 



(A15) 



S„=0 S S =Q 



where the summations run again over both integers and half-odd integers. 

We can perform the summations of Eq. (|A15[) in the integers S and V instead of in S n and S s . Indeed, the first 
factor cancels all the terms with S and T> non-integer so that, 



Nl/2 [JVa/2-S„] 

E E 

S_=0 S s =0 



1 + (-1)2(S„ + S S ) 



Thus wc find, 



■(Srj and S s both either integers or half odd integers) 

Nl/2 +S 

= E] E] ' ' ' an d 2? integers) . 
<s=o T>=-S 



Nl/2 +S 

Mo* - E E ((<$ + l) 2 - 2? 2 ) x 53(5,2?). 



5=0 X>=-S 



The use of the result (|A13[) then leads to, 

Nl/2 s 

Mo t = E E (os+i) 2 -2? 2 

5=0 D=-5 



^v 2 

AT2/2 - 2? 



Nl 

Nl/2 -5 - 1 



^ \ f Nl 
Nl/2 -S) + Wj2 -S-2 

Nl \ ( Nl 
Nl/2-V + l) + \Nl/2-V-l)\ 



Replacing the variable S by 5' = S + 1 we reach a more tractable expression for Aftot, 

Nl/ 2+1 S'-l 

Mot- E E T(S',D), 

5' = 1 U=-5' + l 



where 



T(5', V) = {{S'f - V 2 ) x S(5' - 1, 2?) 

iV 2 



= 5' 2 - 2^ 



Nl/2 - V 
Nl 



N 2 



Nl/2-S' + lJ \Nl/2-S'-l 



N 2 



Nl 



Nl/2-S'J [\Nl/2-V + lJ \Nl/2-V-l 



Nl 



(A16) 



(A17) 



(A18) 
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J E [q(N 2 a +2- q )-p(N! + 2- P )} 



p,q=0 



TV 2 



N 2 

1 a 

p 



p-2 



N 2 a 
p-1 



N 2 

q 



3-2 



This expression can be simplified noticing that, 



TV 
x-2 



Replacing in Eq. (|A20j ) one then finds, 

N^+2 



q{N 2 a +2-q)-p(N 2 a +2-p) 



= -E 

p,q=0 



= -2 



N 2 



N 
x - 1 



iV + 2 

a; 



AT a 2 + 2 
P 



P-1 



TV 2 +2 
9 



p,q=0 

N'?+2 



N 2 



<7=0 

Finally, the use of the identities, 



N 2 + 2 
P 

Nl+2 



p-1 



TV 2 + 2 

q 



p=0 



N~+2 



q=0 



N„+2 



p=0 



N 2 

p-1 



AT 

E 

fe=0 



TV 



= 2 



A' 



(A19) 



is completely symmetric in the summation variables. 

Since T(S', D = ±5') = 0, we can extend the summation over V of Eq. ([AT7f to V = ±S'. We then formally extend 
the summation over S' to S 1 = because the corresponding term vanishes: T(S' = 0,2? = 0) =0. Futhermore, 
T(±5', 2?) = T(S', ±2?) = T(S' , 2?), and due to the symmetry S' <-> 23 we can write, 

E E T(S',V) = l E T(5',23). 

S' = 1V=-S' + 1 S' ,V=-(Nl/2+l) 

Let us introduce the numbers p and g such that, 

S' + iV 2 /2 + l=p & S' =p- (N 2 /2 + 1) 

V + N 2 J2 + l = q ^ V = q- (N 2 /2 + 1) . 
The use of (| A19|) then allows rewriting (|A17[) as, 

Nl+2 



(A20) 



(A21) 



N+2 

E 

fe=0 



TV 
k - 1 



JV+1 

E 

fc=i 



N 
k - 1 



E ( y 



k'=0 



and 



JV 



EK^-fe) 



AT' 



iV-l 

E 

fc=i 



AT! 



(k - l)l(N - k - 1)! 



= N(TV - 1)2 



JV-2 



N~2 



AT (TV — 1) E 



fc-l=0 



TV-2 
fc-1 



N+2 / N \ N+1 ( N \ N f N \ 

E^ + 2-fc)^ ^ = E fc ( iV + 2 - fc )( fe _ 1 j= E ^+2 J 

/c — & — 1 & - 1 — 



JV 



j2(k'+i)(N-k'+i)(l) = e [k'(N-k') + (N +!)](;; 



fc'=0 x 7 fc'=0 

i /at , _ nJV-2 



A^ 



= TV(TV - 1)2^ + (TV + 1)2^ = 2 W " 2 [TV(TV - 1) + 4(TV + 1)] 
= [N 2 + 3N + 4] 2 N - 2 , 
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leads to, 

Mot - \ { [N 2 a + 3N 2 a + 4] 2^- 2 x 2<+ 2 - (N 2 a + 2)(N 2 a + 1)2 N « x 2<} = i2 2N « x 2 = 4^ 2 , (A22) 
which is the desired result. 
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